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ABSTRACT 



The attenuation of surface waves caused by underwater 
turbulence was investigated in a wave-tank experiment. The 
waves studied (frequencies ranging from 1.2 to 12.3 Hertz) 
were strongly attenuated by a zone of grid-generated turbu- 
lence. This attenuation depended on the length of the 
turbulent region, and on the frequency of the incident wave. 
The equation governing attenuation was: 

u u 

h = h e 
o 

where h is wave height, h^ is undisturbed wave height, x is 
the length of the turbulent region, and a is an attenuation 
coefficient proportional to the cube root of frequency. It 
was also noted that the waves were shifted in phase as they 
passed through the turbulence, and that the magnitude of this 
shift increased with frequency. The quantitative results 
of the experiment were obtained from measurements of nearly 
sinusoidal waves. They were confirmed, qualitatively, for a 
continuous spectrum of waves by measurement of a wind-gener- 



ated model sea surface. 
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A LABORATORY EXPERIMENT ON 



SURFACE WAVE ATTENUATION DUE TO 
UNDERWATER TURBULENCE 



1. Introduction 



Objective 

The objective of this experiment was to determine, 
quantitatively, the effect that a zone of underwater turbu- 
lence has on a gravity wave which passes through it. 

P. Groen [l954, p.V^ predicted that a decay in ampli- 
tude and an apparent increase in period of swell would occur 
after a gravity wave train passed through a turbulent region, 
and that this effect is related to external turbulence, that 
is, turbulence which is produced independently of the wave. 
Groen further predicted that this effect depends on both 
wavelength of the waves and characteristic wavelength of 
the turbulence. 

Groen also states that gravity waves may be affected 
by self-produced or internal turbulence. This will not be 
considered in this paper. 

T. P. Barnett [l968j states that one of the possible 
drawbacks of the present scheme of predicting generating 
and dissipation of ocean waves is the neglect of some un- 
known dissipative mechanism. Since his scheme did not 
account for external turbulence, the interactions between 
waves and turbulence may provide the needed explanation for 
loss of wave energy. 
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Because of the difficulties of precise measurement of 
waves and turbulence at sea, this possible energy dissi- 
pation was studied in a laboratory experiment. 

Outline of Experiment 

A wave tank equipped with wave and turbulence generat- 
ing devices was used to produce the desired experimental 
conditions. A traveling surface wave was generated at one 
end of the tank and absorbed by a porous beach at the 
opposite end. The wave passed through a zone of turbulence, 
produced by an oscillating grid in the center of the tank 
bottom. The effect of the turbulent region on the wave 
was evaluated by wave height measurements taken with two 
resistance-type probes located on each side of the turbu- 
lent region. Measurements of turbulent velocities were taken 
with a constant temperature anemometer. 
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2 . Equipment 



Wave Tank 

Dimensions of the rectangular-shaped wave tank used in 
this experiment are shown in Figure 1. 

A false bottom installed in the tank isolated the 
driving mechanism for the oscillating grid and prevented 
it from generating currents which might advect turbulence 
beyond the desired region. The false bottom was stepped 
down in the center of the tank (Figure 1) to accomodate 
the grid which was supported by four metal rods . Since 
the turbulence generator was located entirely below the 
false bottom, no equipment penetrated the water's surface 
except the chain drive for the turbulence generator. A 
false side of the tank prevented this chain from inter- 
fering with the wave as it traveled down the tank. 

The effective depth for the higher-frequency waves 
generated was the depth of the false bottom, approximately 
12 centimeters. When lower frequency waves were generated, 
the false bottom was removed, allowing currents to form. 

The effect of these currents on the experiment appeared to 
be negligible. 

Only deep-water waves were studied to avoid compli- 
cations resulting from the effects of reflection of shallow- 
water waves which would occur in the stepped-down region. 
(Deep-water, in this experiment, refers to waves whose ratio 
of depth to wavelength is greater than or equal to one-half.) 
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Fig, 1 Schematic Diagram of Wave Tank 



Wave Generator 



Waves were generated at one end of the tank by a flat 
aluminum paddle driven by a variable-speed motor. The 
period of the wave generated depended on the speed of the 
motor; the wave height depended on the stroke length of 
the paddle. The equipment used provided two paddle-stroke 
lengths and continuous changes of motor speed. 

Turbulence Generator 

Turbulence was generated by a vertically-oscillating 
grid located in the center of the wave tank just above the 
false bottom. Various screening materials were tried for 
the grid, but it was found that only screening with a 
large, solid surface-to-hole ratio would generate turbu- 
lence strong enough to disturb the near-surface area of the 
water. The grid screen used for this experiment was woven 
of metal strips one centimeter wide, set one centimeter 
apart across the width of the grid and two centimeters 
apart lengthwise. 

The grid was oscillated by apparatus located beneath 
the false bottom of the tank (Figure 2) . All moving parts 
of this apparatus were lubricated by the water which 
surrounded them. A frame, connected to the grid by four 
bolts which penetrated the bottom of the tank, was moved up 
and down by a piston driven through a crankshaft. A chain 
drive connected the crankshaft (through the false side) to 
an electric motor located just above the tank. 
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Fig, 2 Turbulence Generator 
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Wave Probe Design and Calibration 

Resistance-type wave probes were used to measure wave 
height outside the turbulent zone. The probes operate on 
the principle that two parallel wires immersed in water, so 
that the plane of wires is parallel to wave crests, will 
have a resistance that is a function of the depth of immer- 
sion. A block diagram of the probe circuit is shown in 
Figure 3. The power supply is an audio oscillator. The 
probe, which consists of two wires supported by an insulat- 
ing block, is in parallel with one leg of a balanced 
bridge circuit. As the water level varies on the wires, 
the AC signal from the audio oscillator is amplitude mod- 
ulated. The bridge output is amplified, and then converted 
by a demodulator from an amplitude-modulated AC signal to a 
varying DC signal. This signal can then be recorded on 
magnetic tape, strip chart or other suitable recorder. 

Static calibration of the wave probes was accomplished 
by changing the depth of immersion of the probes and plotting 
the water level versus the DC output of the probe circuit. 
Figure 4 shows a typical plot of output versus water level. 
The magnitude of the linear range for the probes used was in 
excess of 7 centimeters. Since the wave heights generated 
were all less than this value, the static response of the 
probes imposed no significant restrictions on the experiment. 

Dynamic calibration of the wave probes was more diffi- 
cult to accomplish than static calibration. The probe and a 
rod, graduated in millimeters, were placed in the wave tank 
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Fig. 3 Wave Probe Circuit 
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Calibration Curve 
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in such a position that water level could be read on the rod. 
When waves were generated in the tank, the difference between 
the maximum and minimum reading on the rod gave the wave 
height. In order to facilitate reading the wave height from 
the rod, a strobe light was used to "stop" the wave train. 

The light flashed at the same frequency as the wave, illum- 
inating only one phase of the wave and thus, eliminating 
the effect of the moving wave on the visual reading. Waves 
of various frequencies were generated and wave height was 
recorded both visually and from the DC calibrated probe 
circuit. The ratio of probe wave height to actual (visual) 
wave height was computed and plotted as a function of fre- 
quency as in Figure 5. The ideal ratio of voltage output 
to wave height was obtained from the static calibration 
and is shown as the value at zero frequency in Figure 5. 

As shown in the figure, height response was excellent for 
frequencies up to 4 Hz and was still fairly good at 26 Hz. 

The absence of readings between 4 and 26 Hz is due to the 
limitations of the wave generators. (Frequencies 4 Hz and 
lower were generated by the afore-mentioned generator; the 
26 Hz wave is a capillary wave of much lower amplitude 
generated by a fan.) 
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Fig. 5- Dynamic Calibration of Wave Probe 



Anemometer 



A calibrated, quartz-coated, hot film probe was used 
with a DISA Constant Temperature Anemometer to measure par- 
ticle velocities in the turbulent water region. A DISA 
Linearizer altered the output of the anemometer so that out- 
put voltage of the circuit was directly proportional to 
velocity of the fluid. 

The anemometer was calibrated by putting the probe 
in a flow of known velocity which was varied and compared 
with output voltage. The values obtained were plotted and 
the curve was corrected for linearity by adjustment of the 
linearizer. A cylindrical tank fitted with an outlet pipe 
near the bottom was used to provide the flow. As the level 
in the tank changes, the outflow velocity varies and a plot 
of velocity versus output is obtained. Values of velocity 
were calculated from the simplified Bernoulli equation: 

U = (2gh)^ 

U is outflow velocity 
g is acceleration due to gravity 
h is height of water in the tank 

A typical calibration curve is shown in Figure 6. 

The use of a constant temperature anemometer in this 
experiment had certain drawbacks: 

The anemometer is normally used to measure small 
velocity fluctuations in a mean flow; the output of the 
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anemometer is the sum of the mean flew and fluctuation: 

U = U + U* 

where U is mean flow velocity, U* is the time varying fluc- 
tuation and U is the flow at any instant. The anemometer 
senses the magnitude of U and thus, in most applications, 
the turbulent fluctuations can be obtained by subtracting 
the mean flow. In this experiment, however, U is zero, or 
very nearly so. Therefore, since the anemometer senses 
only the magnitude of U (U* when U is zero) , as U* changes 
direction, the output of the anemometer remains a positive 
voltage and thus does not give a true velocity profile. 

As a result, there is a tendency for the instrument to 
double the frequency of the fluctuations. 

A second drawback in the use of the anemometer in this 
experiment affects the magnitude of velocities in the tur- 
bulent region. The anemometer probe is most sensitive to 
the component of flow velocity perpendicular to the film. 

When flow in the tank reverses, it passes back over the 
posterior of the probe. The probe, therefore, influences 
the flow physically. When measureing this altered flow, 
the probe is also measuring a flow velocity moving in a 
direction opposite from the flow direction it has been 
calibrated to measure. 

When probe orientation was reversed in a constant flow 
of known unidirectional velocity, probe response was measured 
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at 20 to 25% of the true value for velocities of five to 
six cm/sec. Because of this reduced sensitivity, the probe 
tended to respond primarily to velocities in one direction; 
however, the response in the opposite direction is great 
enough to make interpretation of the probe output very un- 
rertain . 
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3. Spectra of the Paddle-Driven Surface Waves 



The spectral energy of the waves produced by the wave 
generator was strongest at the frequency of oscillation of 
the paddle. However, Figure 7, which is a plot of root- 
mean-square height distribution versus frequency, shows 
that there was also energy at the second and third har- 
monics of the paddle frequency. 

The distribution of wave height with frequency was 
determined as follows: 

1. The wave generator was run until a steady state 
was approximated (i.e., until the waves in the tank no 
longer increased in height) . 

2. Wave height was measured at the two probe positions 
when no turbulence was being generated. 

3. The outputs of the two probes were recorded on a 
PI 6200 magnetic tape recorder operating at .375 inches per 
second . 

4. The tape was then played back at 3.75 inches per 
second and the output was analyzed with a General Radio 
Company Type 1900 Wave Analyzer. (The tape was played back 
at ten times recorded speed so that waves of 1 Hz and greater 
could be measured since the effective low frequency limit of 
the analyzer is 10 Hz.) 

In an application where the signal measured is broad 
band (i.e., for study of the turbulence spectrum) the meter 
of the wave analyzer indicates the value of the instantaneous 
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frequency in Hertz 

Fip-. 7 . RMS Height distribution for Paddle-Driven Waves at Probes #1 & #2 



voltage passed through the filter as averaged by a smooth- 
ing circuit in the equipment. This value must be corrected 
for bandwidth of the filter {the analyzer has a selection 
of 3 , 10 and 50 Hz bandwidths) and converted to an rms 
voltage per one Hertz bandwidth. This voltage was then 
converted to wave height and it is this reading that is 
plotted in Figure 7. 

This method of determining height distribution is 
particularly suitable for presenting broad band signals 
such as noise and turbulence. It does not, however, provide 
an accurate description of the surface wave spectra deter- 
mined in this experiment, since it was found by using band 
widths of 10 and 3 Hertz that the wave energy was essen- 
tially independent of band width and, therefore, the wave 
analyzer voltage reading indicated the value corresponding 
to the true wave height. 

Figure 7 also shows the spectral density of the cir- 
cuit noise. Circuit noise was converted to equivalent wave 
height for comparison with the magnitude of wave energy. In 
most cases, circuit noise was negligible. The peak width of 
the wave spectra in these figures is due to the filter; in 
later figures, the wave spectra will be presented as spikes 
of energy. 

Attenuation due to propagation distance was measured by 
comparing the two probe outputs. In Figures 8 through 14, 
wave height is plotted versus frequency at the two probe 
locations (1.2 meters apart) for the various frequencies 
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Spectral Analysis for a 1.2 Hertz Wave 
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Fig, 9 Spectral Analysis for l.S Hertz 
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Fig* 10 Spedtral Aialysis for a 2,4 Hertz Wave 
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Fig, 12 Spectral Analysis for a 4*1 Hertz 
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Fig. 14 Spectral Analysis for a 4«0 Hertz 



used is this experiment. Figure 15 shows the per cent 
wave attenuation as a function of frequency caused by this 
1.2-meter distance. Figure 15 illustrates that attenuation 
due to distance between the probes was small. In most cases, 
the fundamental frequency of the waves was only slightly 
affected; however, the second and third harmonics were 
significantly reduced in height with distance, that is, 
the wave became more nearly sinusoidal with distance from 
the source. The greater attenuation of the high-frequency 
components is as expected; however, the large amount of 
scatter in the values obtained makes it difficult to assign 
a good relationship between frequency and attenuation in 
this case. 
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4 . Turbulence Spectra 



Typical spectra measured with the constant temperature 
anemometer are shown in Figure 16 for the case of waves, 
waves and turbulence, and turbulence alone. The velocities 
shown were measured at a point in the center of the turbu- 
lent zone, two centimeters below the surface with the probe 
in a position such that the anemometer was most sensitive to 
velocities perpendicular to a vertical plane in line with 
the wave crests. Other similar turbulence spectra were 
measured at six equally-spaced positions at depths of one 
and three cm. 

The drawbacks noted in Section 4 make interpretation 
of the spectra difficult; however. Figure 16 shows signifi- 
cant turbulent energy at all wave frequencies used. The 
curve for waves alone indicated that the probe sensed the 
orbital velocities of the wave and its harmonics, and that 
the high frequency portion of the spectrum of waves and 
turbulence was almost identical to the spectrum of turbu- 
lence alone. The low-frequency portion of the turbulence 
spectrum showed a marked increase in energy at the funda- 
mental frequency of the wave. 
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Frequency in Hertz 

Fig. 16 Velocity Distributions in the Turbulent Region 



5. Attenuation of Waves due to Turbulence 



Attenuation of waves due to turbulence was determined 
by comparing the two wave probe outputs which were recorded 
while turbulence was being generated. Some typical results 
of this comparison are presented in Table 1. Wave energy 
per unit surface area was calculated from the following 
relationship : 

E = 1/8 dgH^ 
d = density 

g = acceleration due to gravity 
H = wave height 

Other comparisons of attenuation due to turbulence 
are shown in Figures 17 to 23. Figure 24 is a plot of 
percent attenuation versus frequency on a log-linear graph. 

The points which are marked by dashed circles represent 
measurements of waves which were of sufficient wavelength 
to undergo a shoaling effect as they passed over the turbu- 
lence generator (in runs with the false bottom removed) . 

The points that fall on the straight vertical line were com- 
pletely attenuated. 

Because of the additional effects which shoaling might 
introduce, the points circled by a dashed line cannot be con- 
sidered good data points for quantitative turbulent attenua- 
tion measurements. They do indicate qualitatively that atten- 
uation is less at lower wave frequencies. 
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TABLE 1 



TYPICAL VALUES OF ENERGY LOSSES 
(Energy density in Ergs/cm^) 



Loss Due to 



Freq . 
(Hz) 


Loss 


Due to Distance 


Turbulence and 


Distance 


Energy : 
Probe #1 


Density 
Probe #2 


% Loss 
of Energy 


Energy 

Probe#l 


Density 
Probe #2 


% Loss 
of Energy 


3.1 


372 


356 


.43 


347 


18.1 


94.8 


* 6.2 


5.63 


1.55 


72.6 


45.4 


.24 


99.6 


** 9.3 


.01 


.00 


— 


.04 


.00 


— 


3.5 


1148 


1127 


.18 


1750 


38.1 


98.0 


* 7.0 


59.3 


34.5 


40.8 


60.2 


.00 


— 


**10.5 


16.5 


1.31 


92.2 


1.59 


.00 


-- 


4.0 


1018 


870 


14.5 


986 


60.00 


93.9 


* 8.0 


22.1 


10.6 


52.0 


21.3 


.00 


— 


**12.0 


.02 


.00 


84.8 


.00 


.00 





* second harmonic 
** third harmonic 
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Fig, 17 Effect of Turbulence on 1.2 Hertz Wave 
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Fig. IS Effect of Turbulence on l.S Hertz 
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Fig. 19 Effect of Turbulence on a 2.4 Hertz Wave 
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Fig, 20 Effect of Turbulence on 3«1 Hertz Wave 
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Fig, 21 Effect of Turbulence on 3.5 Hertz V/ave 
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Fig. 22 Effect of Turbulence on 4.0 Hertz 
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Fig. 23 Effect of Turbulence on 4.1 Hertz Wave 
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Fig* 2^4 Percent Wave Attenuation due to Turbulence 
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The length of the turbulent region used in this experi- 
ment was essentially the same for all data runs. The waves, 
however, were of various wavelengths and thus, their exposure 



to the turbulence differed. A high-frequency wave had more 
wavelength exposed to the turbulence at any one time than a 
low-frequency wave. This unequal exposure to the turbulence 
has been compensated for in Figure 25 where frequency is the 



abscissa and 



Ah/h 



Ax/L 

length, is the ordinate. 



, the fractional change of height per wave- 



where; H = h-h^^ 

h = height of the incident wave corrected for 
attenuation due to distance by multiplying 
the value obtained at the first probe by 
(1-b) where b is the percent attenuation due 
to distance alone for that frequency, 
h^ = height of the wave after passing through the 
turbulent region. 

Ax = length of the turbulent zone, and 
L = wavelength 

Figure 25 suggests that the relationship between log 
and log (frequency) is linear. Assuming this to be so, a 
least squares fit to the points results in a line given 
approximately by the equation: 

log ^ 
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Fig. 25 Attenuation of Waves Due to 
Turbulence Plotted as a Function of 
Number of Wavelengths in the Turbu- 
lent Region 
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which can be written 



< ^h _ Kf-— 



( 2 ) 






where K has units of time to the^f ive-thirds power and a 
numerical value of approximately 1.26. 

If the waves are assumed to be in deep water, then: 



L = 



21Tf‘ 



and equation (2) can be written as 



h - q 2lTKf 



(3) 



If Ah and AX are allowed to become very small, equation (3) 
can be written in differential form: 



^ . dx fV3 

h q 



(4) 



Integrating both sides of equation (4) yields: 



h = h e 
o 



-ax 



(5) 



where h is the initial height of the wave and a is an 

2'irK .q 

attenuation coefficient which is equal to — ! — f t ^nd has 

g 

units of reciprocal length. 

This relation may be expressed in terms of energy 
density as follows: 

E = |- dgh^ 



where d is fluid density and E is energy per unit area. 
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or 



where 



Substituting equation (5) : 



_ 1 j v,2 <=“2ax 

E = g dgh- e 



-2ax 



E = E e 
o 



E = - dgh ^ 
o 8 o 



(X, ' 



‘iT 



K ^ 



'b 






5 



(6) 



'/!( 
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6. Phase Shift in the Turbulent Zone 



The fundamental frequency of the wave trains experienced 
a phase shift upon passing through the turbulent region. 

This shift was not measured for second and third harmonics 
because of the limitations in the measuring technique. 

Phase shift was measured by recording the output of both 
probes on a two-channel strip chart recorder. Waves alone 
were measured initially; the turbulence generator was then 
turned on and the outputs of the two probes were compared 
with respect to their phase relationship both before and 
after passing through turbulence. 

Firure 26 is a plot of frequency versus phase shift and 
frequency versus phase shift per number of wavelengths pass- 
ing through the turbulent region at any one time. 

Figure 26 shows that at the lowest frequency measured, 
the phase shift was negative. At higher frequencies the 
phase shift was positive and increased with frequency. 

Assuming that the frequency of the waves remained con- 
stant as they travelled the length of the tank, the observed 
phase shift implies that both the wavelength and the phase 
speed of the waves changed as they passed through the turbu- 
lent region and that this phase change increased in magni- 
tude with frequency. It was observed visually, particularly 
for the higher frequency waves , that as the waves entered the 
turbulent region, they steepened considerably at first and 
were then reduced in height and steepness in the 
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Q Phase shift in radians/no, of wavelengths passing 
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Fig* 26 Phase Shift Versus Frequency 



later portion of the turbulent region. 



It is uncertain in 



what manner, 
to the phase 



if at all, this observed steepening is related 
shift . 
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7. Attenuation of Wind-Driven Waves 



A model sea surface was generated with a small centri- 
fugal blower, and measurements of surface spectra were taken 
with and without the turbulence. Figure 27, a plot of rms 
wave height versus frequency, shows the results of these 
measurements. The entire water surface was in the blower's 
generating area. Thus, the waves were not only affected 
by propagation distance and turbulence, but also by being in 
a region of generation. The figure indicates that the higher 
frequency waves were almost completely attenuated, while a 
small amount of energy remained in the lower frequencies. 

The difference in peak values between the undisturbed spectra 
at probes one and two shows the effect of a longer fetch for 
the second probe. The blower also generated many high fre- 
quencies which were attenuated by propagation distance before 
they reached the second probe. 
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Fig. 27 Attenuation of the Model Sea Surface 



8 . Summary and Conclusions 



This experiment showed that underwater turbulence had 
the following effects on surface waves: 

1) Wave height decreased as a result of attenuation 
by turbulence for all frequencies measured (1.2 to 
12.3 Hz) according to the formula 

v^ _ v^ 

h = h e 
o 

-1 / 

where a is an attenuation coefficient equal to — — f 

g 

2) The waves studied were shifted in phase by turbu- 
lence; below 1.5 Hertz, this shift was negative and 
became positive above 1.5 Hertz. 

3) Measurements using a model sea surface indicate, 
qualitatively, that a continuous spectrum of waves 
experiences attenuation in much the same manner as 
the nearly sinusoidal waves of this experiment. 
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